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The evolution of weak discontinuity is investigated in the flat FRW universe with a single scalar 
field and with multiple scalar fields. We consider both massless scalar fields and scalar fields with 
exponential potentials. Then we find that a new type of instability, i.e. kink instability develops in 
the flat FRW universe with massless scalar fields. The kink instability develops with scalar fields 
with considerably steep exponential potentials, while less steep exponential potentials do not suffer 
from kink instability. In particular, assisted inflation with multiple scalar fields does not suffer from 
kink instability. The stability of general spherically symmetric self-similar solutions is also discussed. 
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I. INTRODUCTION 

Future precision observations of the cosmic microwave 
background radiation and/or gravitational lensing phe- 
nomena will provide a large number of new data to be 
analyzed and interpreted. The standard cosmology as- 
sumes that the observable part of the universe on scales 
of the order of the present Hubble distance is almost spa- 
tially homogeneous and isotropic. Actually, the data of 
WMAP for the cosmic microwave background radiation 
tell us that the universe at the last scattering surface is 
well approximated by the Friedmann-Robertson- Walker 
(FRW) universe so that understanding of the detailed 
properties of the FRW universe is an important step in 
cosmology 

One of the most important problems is how cosmo- 
logical structures on various scales formed from the very 
homogeneous and isotropic universe at early times. It 
is known that the gravitational instability scenario with 
cold dark matter well describes structure formation from 
small primordial perturbations in our universe 

Recently, a new kind of instability has been found in 
the studies of self-similar solutions. Self-similar solutions 
are defined in terms of a homothctic Killing vector field in 
general relativity and include a class of the flat FRW 
solutions. The instability concerns weak discontinuity. 
If the perturbations with weak discontinuity are inserted 
into a self-similar solution which is unstable against this 
mode, the discontinuity grows as time proceeds. This in- 
stability has been called kink instability. Kink instability 
of self-similar solutions was originally found by Ori and 
Piran in spherically symmetric isothermal gas systems in 
Newtonian gravity [^|. In general relativity, Harada in- 
vestigated the kink instability of self-similar solutions for 
the spherical system of a perfect fluid with the equation 



of state p = k/i [3- The present authors investigated the 
kink instability of self-similar solutions for the spherical 
system of a stiff (k = 1) fluid and those of a massless 
scalar field 6]. These works have shown that the kink 
instability can occur in a large class of self-similar solu- 
tions. In particular, the work in general relativity showed 
that the flat FRW universe is unstable for a perfect fluid 
for 1/3 < k < 1 and for a massless scalar field. 

In this letter we derive a stability criterion against kink 
mode perturbations for the power-law flat FRW universe 
with a single scalar field and with multiple scalar fields 
with and without exponential potentials, which arise nat- 
urally in supergravity Q or theories obtained through 
dimensional reduction to effective four dimensional theo- 
ries 0, 0] . The stability of general spherically symmetric 
self-similar solutions with scalar fields is also discussed in 
Appendix lAl We adopt the units such that c = G = 1 
and the abstract index notation of [Icj . 



II. STABILITY CRITERION FOR THE FLAT 
FRW UNIVERSE 



power-law flat FRW solution with and without 
exponential potentials 



We begin with the action which describes the self- 
gravitating system of scalar fields with and without ex- 
ponential potentials. 



S = 
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(2.1) 

where {Vi\i — 1, •■• ,n} and {Xi\i — 1, ■■• ,n} are n 
real non- negative and positive constants, respectively. 
In this model we currently have n real scalar fields 
{4>j\j = 1, • • • , n}, for which the total potential has the 
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form 



where {/%} are constants satisfying 



Vtc 



(2.2) 



If Vj = 0, then <j>j is massless and Xj is meaningless. 
We assume < Ai < A2 < • • • < A n without loss of 
generality. This action has been considered in the context 
of assisted inflation 0, 0, . The action involving a 
perfect fluid has been studied by several authors 

Then the total energy-momentum tensor is given by 



T ab = (V^VVi - | 5 a "V c ^V c ^) - g ab V t 



i=l 



while equations of motion for {<f>j} is given by 

dV tot 



(2.4) 



We consider the power-law flat FRW spacetime with 
the metric; 



ds 2 



~dt l 



(dr 2 +r 2 dfl 2 ), (2.5) 



where a and to are constants and dti 2 = dQ 2 + sin 2 Qdip 2 . 

We first show that there are no power-law flat FRW 
solutions which include both nontrivial massless scalar 
fields and scalar fields with exponential potentials. When 
there exist m scalar fields {<^-|l < j < m} with a total 

potential J^fcLi ^ ex P( — VSTrXi4>i) and (n — m) massless 
scalar fields {4>j\m + 1 < j < r?,}, the Einstein equations 
Gab = 8TrT a b and the equations of motion for scalar fields 
are written as 



3a(a - 1) = 87rf 2 



E' 



a(3a — 1) 



m 

8nt 2 Vi 



VWirXi 



/8TT\i<t>i 



(2.6) 
(2.7) 



i=l 



t 2 <j>j,tt 



3at(f> 



j,t 



y/&K\ j V j t 2 e-^ X i't>i 

(for 1 < j < m),{2. 
(for m + 1 < j < n), 



where the comma denotes the partial derivative. From 
Eqs. l|2~^|l and |2"T)l . we obtain 



a = 47r<2 E^M' 

i=\ 

so that {4>j} must have the form 

7 - = — k, lni + C; 



(2.9) 



(2.10) 



a = 4tt^k 2 , 



(2.11) 



and {Cj} are constants. We can set Cj ~ for all j 
by redefining the scalar field. From Eq. (|2 . T|> and the 
assumption of Vj > for 1 < j < m, Kj = — 2/(v / 87t\\j) 
must be satisfied for 1 < j < m. Then Eqs. I|2.7|) and 
(|2.8|) reduce to 



a 



2 

Et2' 



(2.12) 



(2.3) On the other hand, we obtain from Eq. I|2.11|l that 



r-1 A? 

4 — 1 1 



i—m-\-l 



(2.13) 



Eqs. I|2.12|l and l|2.13|l give a contradiction Kj = for m+ 
1 < J < ^i- Therefore, we conclude that if the power-law 
flat FRW solution includes both massless scalar fields and 
scalar fields with exponential potentials, those massless 
scalar fields must be constant, i.e., trivial. 

When there exist only massless scalar fields, the Ein- 
stein equations and the equations of motion for scalar 
fields give a — 1/3 and 47r^™ =1 K 2 ~ 1/3, while when 
there are no massless scalar fields, they give 



E \2' 



Vj 



(2.14) 
(2.15) 



The latter equation gives a condition ^2™ =1 \~ 2 > 1/6 
for Vj > 0. Setting n = 1 in both cases, we can obtain 
the results for a single scalar-field case. 

For the later convenience, we transform the metric 
()2.5|l to the self-similar form l|B2jl . As shown in Ap- 
pendix El the flat FRW solution is self-similar if and 
only if the scale factor obeys a power-law. Therefore 
it is concluded that there are no self-similar flat FRW 
solutions if there exist both nontrivial massless scalar 
fields and those with exponential potentials. Because 
the case with a = 1 is exceptional, in which the homo- 
thetic Killing vector is not tilted [l5j], we consider only 
the case with a / 1. With a coordinate transformation 
r = 1 1 — al^tpf 1- ", the metric (|2.5|l is transformed to 
the self-similar form (|_B2|) as 



ds 2 = -dt 2 + z 2a (df 2 + (1 - a 



- 2 r 2 dn 2 ) 



(2.16) 



where z = t/f is the self-similarity coordinate. 

The particle horizon z = 1 plays an important role in 
the following analysis. This is because it is a similarity 
horizon, which is given by a radial null ray on which z is 
constant. 
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B. kink instability 

We consider general full-order spherically symmetric 
perturbations with scalar fields such as 

ds 2 = _ e 29(t,f) dt 2 + z 2 ae 2u,(t,f) df 2 

+ {l-a)- 2 z 2a f 2 e 2n{t ' f) dQ 2 , (2.17) 
4>j = -Kjhit + ipj(t,f). (2.18) 

Hereafter we adopt the isotropic coordinates, in which 
Q — fj. We define the following coordinates: 

Z = hxz, (2.19) 
T = hit, (2.20) 

in which Eqs. (|2.17() and i|2.18[l are represented by 

ds 2 = „( e 2-(T,Z) _ e 2( Q -l)Z + 2^(T,Z) )e 2T dT 2 

_ 2e 2( a -l)Z+2T+2u, dTdz + e 2( a -l)Z+2T+2u> dz 2 

+ (1 - a)- 2 e 2{a -^ z+2T+2 ^dVL 2 , (2.21) 
<P d = ~ Kj T + ^ 3 (T,Z), (2.22) 

where 

a(T,Z) = a(t(T),f(T,Z)), 

uj(T,Z) = Cu(t(T),f(T,Z)), 

i>.j(T,Z) = ^{t{T),f(T,Z)), 

for 1 < j < n. 

We consider perturbations a, uj and {ipj} which satisfy 
the following conditions: 

(1) The initial perturbations vanish outside the similar- 
ity horizon. 

(2) a, lo, {ipj} and {tp'j} are continuous, in particular at 
the similarity horizon Z = 0, where a prime denotes the 
derivative with respect to Z. 

(3) {ipj} and {ip'J} are discontinuous at the similarity 
horizon, although they have finite one-sided limit values 
as Z — > —0 and Z — > +0, where a dot denotes the deriva- 
tive with respect to T. 

(4) The quasi-local mass is continuous, i.e., there are no 
singular hypersurfaces in the spacetime. 

Now we consider the behavior of the perturbations 
at the similarity horizon. The perturbations satisfy 
ipj = ipj = for all j and ip'J ^ for some j at the 
similarity horizon at the initial moment T = Tq due to 
condition (2). The evolution of the initially unperturbed 
region is completely described by the background flat 
FRW solution because information from the perturbed 
side cannot penetrate the unperturbed side due to con- 
dition (1). Then, we find a = to = 0, ipj = ip'j = for 
all j and tp'J ^ for some j at the similarity horizon for 
T > T due to conditions (2) and (3). The Misner-Sharp 
quasi-local mass m is defined by 

m = ^(l-R,^), (2.23) 



where R is the circumferential radius. Because of 

R = |1 -ar 1 exp((a- 1)Z + T + uj) 

in the present case, to' = is satisfied at the similar- 
ity horizon due to condition (4). Then we find a' = 
from the equation of motion for scalar fields. From 
the (00)+(01) and (11)+(01) components of the Einstein 
equations 

Rab = 8tt \ T ab - ^9abT^j 

, ui" = and a" = are obtained, respectively. Dif- 
ferentiating the equation of motion for scalar fields with 
respect to Z and estimating both sides at the similarity 
horizon, we obtain 

i,"j - (1 - 2a) V; = \ hm (l - e 2 ^ z ) ^ , (2.24) 

for each j. We can show that the right-hand side van- 
ishes 0] . Finally, the full-order perturbation equation for 
ipj at the similarity horizon is obtained as 

$' = (l-2a)V£. (2-25) 

It should be noted that the perturbations are those of 
full-order although this equation is linear. This equation 
can be integrated to obtain 

ip'j oc e^~ 2 ^ T . (2.26) 

Therefore, it is found that the perturbation decays expo- 
nentially for a > 1/2, it is constant for a = 1/2 and it 
grows exponentially for a < 1/2. 

It is noted that these perturbations are gauge- 
independent as shown in Appendix El Here we define 
instability by the exponential growth of discontinuity. 
Then we find the following criterion: the flat FRW uni- 
verse with a > 1/2 are stable against the kink mode, 
while those with a < 1/2 are unstable. Solutions with 
a = 1/2 are marginally stable against this mode. 

Although we concentrate on the flat FRW solutions 
in this letter, we can obtain the stability criterion for 
general spherically symmetric self-similar solutions with 
regular similarity horizons adopting the similar analysis 
as that in 0, which is summarized in Appendix lAl 

III. DISCUSSIONS 

We have derived the stability criterion for spherically 
symmetric self-similar solutions with scalar fields with 
independent exponential potentials to suffer from kink 
instability. It can be applied to any spherically symmet- 
ric self-similar solutions with regular similarity horizons. 
The kink instability, which we have considered here, was 
studied in general relativity for a perfect fluid with an 
equation of state p = k\x an d for a massless scalar 
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field @. We focus on the flat FRW universe in this letter, 
for which result is summarized in Tables HI and ITD This is 
a new kind of instability of the FRW universe. We first 
discuss the relations of our result to other studies in the 
literature on the system with scalar fields with or without 
exponential potentials. After that, we discuss astrophys- 
ical or cosmological implications of the kink instability of 
the flat FRW universe. 

Let us see the flat FRW universe filled with a single 
matter field. The FRW universe is unstable for a per- 
fect fluid with 1/3 < k < 1, for a massless scalar field 
and for a scalar field with an exponential potential for 
4 < A 2 < 6. These results suggest that if there exists 
a phase in the flat FRW universe in which a scalar field 
dominates other matter fields, the scalar field with an 
exponential potential for < A 2 < 4 is preferred rather 
than a massless scalar field. Such a potential naturally 
arises in supergravity Q or theories obtained through 
dimensional reduction to effective four dimensional the- 
ories nn. 

Kitada and Maeda 0, 0] showed the cosmic no- 
hair theorem, which is a generalization of the result by 
Wald and states that all ever-expanding homoge- 
neous models which contain a scalar field with an expo- 
nential potential together with a matter field satisfying 
the dominant energy condition asymptote to the infla- 
tionary flat FRW solution for < A 2 < 2. According to 
the criterion obtained in this letter, the flat FRW solution 
with 4 < A 2 < 6 suffers from kink mode perturbations, 
which are kinds of inhomogeneous perturbations, so that 
our result does not affect the cosmic no-hair theorem. 

Next let us move on to the case where the flat FRW 
universe is filled with multiple scalar fields. In such 
a situation, Liddle, Mazumdar and Schunck proposed 
assisted inflation, where an arbitrary number of scalar 
fields with independent exponential potentials evolve to 
the inflationary flat FRW solution even if each indi- 
vidual potential is too steep to support inflation by 
its own 

'H3 • Coley and van den Hoogen showed that 
the assisted inflationary solution is a global attractor 
in the spatially homogeneous and isotropic universe if 
ELiV 2 > V 2 Together with the result on the 
spatially homogeneous scalar-field cosmological models 
by Billyard et al. 0], they concluded that the assisted 
inflationary solution is a global attractor for all ever- 
expanding spatially homogeneous cosmological models 
with multiple scalar fields with exponential potentials 
provided Yl7=i A^ 2 > 1/2 13]. According to the cri- 
terion obtained in this letter, the expanding flat FRW 
solution with Y^i=i -V > 1/4 does not suffer from kink 
instability, so that our result does not affect the generic- 
ity of assisted inflation even if we allow initial data with 
weak discontinuity in the second-order derivative of the 
scalar field. 

For any case of the flat FRW universes we have seen 
in this letter, we find that the stability of the flat FRW 
universe against kink mode perturbation changes from 
stability to instability as the power index a of the scale 



factor increases through 1/2, which is clear from Tables [i] 
andllTI This suggests that the cosmological kink instabil- 
ity is related to the expansion-law of the universe rather 
than the characteristics of the matter fields. 

The kink instability is particularly important for the 
flat FRW universe in the radiation (k = 1/3) dominated 
era. It is critical for some cosmological models which 
assume a stage when the effective equation of state is 
"harder" than the radiation fluid. Because the flat FRW 
solution is homogeneous, we can choose any point as a 
symmetric centre in the present analysis. This implies 
that a bubble-like structure will develop and the bub- 
ble walls collide with each other. This consideration 
suggests that density perturbations can arise in differ- 
ent scales as a result of multiple bubble wall collisions. 
Although the product of kink instability has not been 
clear yet, it is plausible that it will bring shock- wave for- 
mation, which could be important for the cosmological 
structure formation. Actually, the shock-wave formation 
has been observed in the numerical simulations of the 
gravitational collapse of the density fluctuations in the 
radiation-dominated flat FRW universe |2p|. The out- 
come of the instability and the effect on the structure 
formation should be comprehensively investigated. 
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APPENDIX A: STABILITY CRITERION FOR 
GENERAL SELF-SIMILAR SOLUTIONS 

Here we show that kink instability may occur in a very 
large class of spherically symmetric self-similar space- 
times. The analysis for the case of a single massless scalar 
field is straightforwardly extended to the case of mul- 
tiple scalar fields. 

We adopt the Bondi coordinates for spherically sym- 
metric spacetimes as 

ds 2 = -ggdu 2 - 2gdudR + R 2 dn 2 , (Al) 

where g = g(u, R) and g = g(u, R). For later convenience 
we define new functions {hj(u, R)} as 

4>j(u,R) = hj(u,R) - Kj]n\u\, (A2) 

where {Kj\j = l,--- ,n} are constants. We define the 
following self-similar coordinates: 

x= , X = ln\x\, U = -\n\u\. (A3) 
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We refer to u < and u > as early time and late time, 
respectively. As u is increased, U increases in early times, 
while U decreases in late times. As u is fixed and R is 
increased, X increases. 

Then the Einstein equation and the equations of mo- 
tion for scalar fields reduce to the following partial dif- 
ferential equations: 



(ln<?)' 




(A4) 



-9, 



(A5) 







(f) 









Tl 



+2x 



{x(hi + % + Ki) -gh'ij, (A6) 
8ttAj gr 2 \u\ ^Vje~ ~ h ' 
(hj + h'j) + (fij + h'jY + Kj] , (A7) 



where the dot and prime denote the partial derivatives 
with respect to U and X, respectively. 



1. self-similar multi-scalar- field solutions 

For self-similar solutions, we assume that g — g(x), 
g — g(x) and hj — hj(x) for all j. From Eq. I|A7(1 . Kj is 
arbitrary if <ftj is massless [21J, while Kj = — 2/(v / 87rA. ) ) 
should be satisfied otherwise for consistency [22j. Drop- 
ping the terms with dots in Eqs. (|A4I) - (|A7() . we obtain a 
set of ordinary differential equations for g, g and {hj}. 

These equations are singular at g = 2x, which is called 
a similarity horizon. We denote the value of x at the 
similarity horizon as X( s ) and also the value of X as 



X (s) = In | 
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We consider self-similar solutions with 



finite values of functions g, g, {hj} and {hj} and their 
gradients with respect to X at the similarity horizon. 
Then we find at the similarity horizon X = Xr s \ 



5(s) 



9(s) 



'•(s) 



(2 + 87T ££=!«?) 



1 

= 2.x 



^xf s) EILi Vie 



( B ), 



7, _ 2K J 

11 m - — 



87rA 3 x (s)g(s) ^-e-^^^M 



(A8) 
(A9) 
,(A10) 



for Kj ^ 0, where the subscript (s) denotes the value at 
the similarity horizon. When Kj — for all j, which can 
be satisfied in the case that all scalar fields are massless, 
i.e. Vj = for all j, we find 



at the similarity horizon. A similarity horizon corre- 
sponds to a radial null curve because 



dX 
UJ 



1 



9_ 

2x 



(A12) 



is satisfied along a radial null curve. No informa- 
tion propagate inwardly beyond the similarity horizon 
in early-time (late-time) solutions for g( s \g~t$\ > (<)0, 
while no information propagate outwardly beyond the 
similarity horizon in late-time (early-time) solutions for 
9(s)9(s) > (<)0. 

A class of the expanding flat FRW solutions belongs 
to the late-time self-similar solutions with an analytic 
similarity horizon. The flat FRW solutions are char- 
acterized by the parameters {Kj}, which are gauge- 
independent. Actually, the FRW solution in the same 
coordinates with those used in this appendix was ob- 
tained by Christodoulou for the case of a massless scalar 
field M. 



2. kink instability of self-similar solutions 

We consider perturbations which satisfy the following 
conditions in the background self-similar solution: 

(1) The initial perturbations vanish inside the simi- 
larity horizon for early-time (late-time) solutions with 
S(s)<7(s) > (<)0- (Conversely, the initial perturbations 
vanish outside the similarity horizon for late-time (early- 
time) solutions with <7( s )<7(s) > (<)0.) 

(2) g, g, {hj} and {hj} are continuous everywhere, in 
particular at the similarity horizon. 

(3) {hj} and {h"} are discontinuous at the similarity 
horizon, although they have finite one-sided limit values 



as X -> X 



(s) 



Oandl^l 



(s) 



0. 



We denote the full-order perturbations as 

6g(U,X) = g(U,X)-g (b) (X), 
6g(U,X) ee g(U,X)-g (b) (X), 
Sh.iU^X) ee hj{U,X)-h j{h) {X), 



(A13) 
(A14) 
(A15) 



where 9(b)? 9(b) an d {^'(b)} denote the background self- 
similar solution. 

Now we consider the behavior of the perturbations at 
the similarity horizon, of which we mean X — > X^ — for 
9(s)9(s) > 0, while X -> A" (s) +0 for 5(s)5(s) < 0. Applying 
the similar analysis as that in [6(, we finally obtain the 
full-order perturbation equation for Shj at the similarity 
horizon is obtained as 



Shj = I -I + Stt^k? ] Shj. 



This equation can be integrated to obtain 



(A16) 



9( S ) = 3(s) 



2x 



O0> h 'i{s) 



h', 



i(s) ' 



(All) 



Shj oc e 



0U 



(A17) 
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where 

n 

13= -1 + 8tt^k 2 . (A18) 

i=l 

Therefore, for early-time solutions, it is found that the 
perturbation decays exponentially for 87r^™ =1 nf < 1, it 
is constant for 8ir Yn=i K ? = 1 an d ^ grows exponentially 
for 87T^™ =1 «| > 1. The situation is reversed for late- 
time solutions. 

It is noted that these perturbations are gauge- 
independent as shown in Appendix [U] Here we define 
instability by the exponential growth of discontinuity. 
Then we find the following criterion: for early-time so- 
lutions, the solutions with regular similarity horizon and 
87r^™ =1 n 2 < 1 are stable against the kink mode, while 
those with 87r Yl7=i K f 1 are unstable. Solutions with 
87r ^™ =1 n 2 = 1 are marginally stable against this mode. 
The situation is reversed for late-time solutions. Setting 
n = 1, we reproduce the result obtained in |6| for the 
case of a massless scalar field. 



These equations result in 

(!)X§g)'© 2 <b*> 

This equation can be integrated to obtain 

far 5 , (B9) 

where s = C(l)/S(l). Then Eq. JB7|) gives 

a 2 (f) oc f 2(1 - 1/s) , (BIO) 

and therefore we conclude that the scale factor obeys the 
power-law: 

a(t) cxi (1 ~ 1/s) . (Bll) 

Inversely, as shown in § III Al we prove that if the 
flat FRW spacetime has a scale factor which follows the 
power-law, then the spacetime is self-similar with a 1. 

It is noted that the case with a = 1 is also homothetic, 
where the homothetic vector is given by td/dt |lfi| . 



APPENDIX B: POWER-LAW FLAT FRW 
SPACETIME AS A SELF-SIMILAR SPACETIME 

We prove that if the flat FRW spacetime is self-similar, 
then the scale factor follows a power-law. We show that 
the flat FRW spacetime with a metric 



ds A = -dt l + a(t) 2 (dr 2 + r 2 dn 2 ), 



(Bl) 



can be transformed to a general spherically symmetric 
metric 

ds 2 = -A(z) 2 dt 2 + B(zfdf 2 + C(z) 2 f 2 dn 2 , (B2) 
with a homothetic Killing vector 

(B3) 



t— -— 

dt df : 



where z = t/f, if and only if the scale factor a(t) obeys 
the power-law. We consider coordinate transformations 
from Eq. i|Bl|) to a self-similar form (|B2|I . t remains 
unchanged, apparently, while we define a new coordinate 
f = f(r) to obey 



a(t) 2 dr 2 = B(z) 2 df 2 



and 



a{t) 2 r 2 = C{zff 



(B4) 



(B5) 



Using a relation t = zr and setting z = 1, we can rewrite 
these equations as 



and 



a 2 (f) = B{l) 2 (df/dr) 2 , 
a 2 {¥) = C(l) 2 (f/r) 2 . 



(B6) 
(B7) 



APPENDIX C: GAUGE-INVARIANCE OF KINK 
INSTABILITY 

To construct the gauge- invariant quantity we introduce 
the 2+2 split of spherically symmetric spacetimes accord- 
ing to [ill ■ We write the spherically symmetric spacetime 
as a product manifold M. = M 2 x S 2 with metric 



= diag(su B ,r 2 7 6), 



(CI) 



where gAB is an arbitrary Lorentz metric on M 2 , r is 
a scalar on M 2 with r — defining the boundary of 
M , and j a i, is the unit curvature metric on S 2 . We 
introduce the covariant derivatives on spacetime M, the 
subspacetime M 2 and the unit sphere S 2 with 



g^u-\ = 0, 
9ab\c = 0, 
gab-.c = 0. 



(C2) 
(C3) 
(C4) 



We also introduce totally antisymmetric covariant ten- 
sors on A4, M 2 and S 2 , respectively: 

£fj,u — ^ffij ^ixv;X = 0, £/j,\£ = 3^, (C5) 



£ab — —£ba, £ab\c — 0, £ac£ 

Sab 



BC 



c c He- cr^ c — rfi 

&baj ^ab:c c ac c g a . 



(C7) 



The kink mode concerns the density gradient on the 
comoving slicing for the fluid case, which is gauge- 
invariant [3, |q|. For the scalar field case, it concerns 
the second-order derivative of the scalar field 0. From 
the second-order derivative, we can construct a similar 
quantity / such as 



I = e 



AB „CD , 



V \B<Pj\D(Pj\AC, 



(C8) 
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for some j. Since the metric, the scalar field and their 
first derivatives with respect to T and X are not per- 
turbed at the similarity horizon for the kink perturba- 
tion, the full-order perturbation 51 = 1 — /(b) becomes 



Therefore, the quantity / has a nontrivial value /(b) at 
the similarity horizon, i.e. 



/ ( b) =e 2U f(X), 



(C15) 



51 



,AB CD 



g <t>j\B<i>j\DHi\ 



AC, 



(C9) 



where 5<pj = <f>j — 4>j(h) ■ 

For the perturbed FRW universe given in § [nj the line 
element on M 2 is 



ds 2 



M- 



_^ e 2a(T,Z) _ e 2(a-l)Z+2u(T,Z)y2T d j,2 
_ 2e 2( a -l)Z+2T+2^ dTdz 

_ e 2{ a -l)Z+2T+2^ dz 2^ (C1Q) 



Therefore, it is clear that for the background flat FRW 
universe, the quantity / has a nontrivial value /(b) even 
at the similarity horizon, i.e. 



(Cll) 



SI in the limit to the similarity horizon is calculated to 
be 



r t 2 -2T i II 

SI = Kj e ^. (s) . 



(C12) 



The factor e~ 2T is the same as that for the background 
flat FRW universe. Therefore, we obtain for the pertur- 
bation evolution in the limit to the similarity horizon 



SI 
h 



(C13) 



For more general case in Appendix the line element 
on M 2 can be written as 



j 2 -IV 

as M 2 = gxe 



2 - - ) dU 2 + 2dUdX 



(C14) 



where f(X) is a function of X. 

In terms of the coordinates (U,X) presented in Ap- 
pendix ^ SI in the limit to the similarity horizon is cal- 
culated to be 



a 2U 



51 = -- 



(C16) 



Or, using the notation presented in Appendix lAl we ob- 
tain 



51= 



(C17) 



in the limit to the similarity horizon. The factor e 2U is 
the same as that for the background solution. Therefore, 
we obtain for the perturbation evolution in the limit to 
the similarity horizon 



f «^ (i) (E0. 



(C18) 



In the whole spacetime M. the quantity / can be writ- 
ten as 

/ i "'// A ^' /;J ,o, ; , ! o, :(iA . (C19) 
It is clear that this is gauge-independent. 
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TABLE I: Stability of kink mode perturbation in the flat 
FRW universe with a perfect fluid with the equation of state 
p — k/i. The power index a of the scale factor is given by 



+ k)]. See g 


for details. 


Parameter 


Power index a 


Stability 


< k < 1/3 

k = 1/3 
1/3 < k < 1 


1/2 <a< 2/3 

a = 1/2 
1/3 < a < 1/2 


Stable 
Unstable" 
Unstable 



"In this case, although the linear stability is marginal, nonlinear 
effect causes nonexponential instability. 



TABLE II: Stability of kink mode perturbation in the flat 
FRW universe with scalar fields. For exponential poten- 
tials, we consider the sum of independent potentials, V to t = 
J2i Vi exp(— \/8n\itf)i), for scalar fields {<f>i}, where the power 
index a of the scale factor is given by a — 2 ^ 4 A" 2 ! . 



Potential 


Parameters 


Power index a 


Stability 


Massless 




a = 1/3 


Unstable 


Exp 


E 4 at 3 = i/4 

l/6<EA- 2 <l/4 


a > 1/2 
a = 1/2 
1/3 < a < 1/2 


Stable 
Marginal 
Unstable 



